Electromagnetic Field Quantization in Time-Dependent Dielectric Media 
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We present a Gupta-Bleuler quantization scheme for the electromagnetic field in time-dependent dielectric 
media. Starting from the Maxwell equations, a generalization of the Lorentz gauge condition adapted to time 
varying dielectrics is derived. Using this gauge, a Gupta-Bleuler approach to quantize all polarizations of the 
radiation field and the corresponding constraint condition are introduced. This new approach is difi'erent from 
the quantized electromagnetic field in vacuum in the sense that here the contributions of unphysical photons 
cannot be thoroughly eliminated, which further lead to a surface charge density. Finally, a discussion of potential 
experimental tests and possible implication is also made. 
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I. INTRODUCTION 

In quantum field theory, attempts to quantize electromag- 
netic field in material media date back to more than half a 
century ago fll, Qj- Recently, interests in thisproblem are 
mainly related to the flourishing cavity QED H and exper- 
iments of quantum optics in media |5|l- In order to solve 
this problem, various methods have been proposed. Of earlier 
works, the quantization scheme given by Jauch and Watson in 
1948 is taking into account only nondispersive and uniform 
media |2]. This scheme has been applied to inhomogeneous 
nondispersive cases too [5]. In recent works, quantization of 
the electromagnetic field in dispersive media, which can be 
either homogeneous or inhomogeneous, has been considered 
In particular, there has been attempt to study the most 
general case with nonlinear materials [9]. More lately, this 
study has also been extended to anisotropic media lUOll . But 
few of these previous works have taken into account the time- 
varying properties of the medium. 

During the last two decades, it was found that a time- 
varying dielectric permittivity may lead to generate quanta of 
electromagnetic field (photons) even from vacuum lll-l-^. 
This phenomenon is similar to pure quantum effects such as 
dynamical Casimir effect (DCE) 1,1 5,1 . For this purpose, sev- 
eral experimental schemes have been suggested by the MIR 
proposal 1 16]. In one of these schemes, one boundary of a 
cavity is made of time-dependent media. So, to clarify how 
to quantize the electromagnetic field in time-dependent media 
will help to understand these experiments. For this purpose, 
many literatures focused on one-dimensional case [12]. They 
split up directly the electromagnetic field into two indepen- 
dent scalar fields assuming to the decoupling of two polariza- 
tions 1 17]. However, this method cannot be directly extrapo- 
lated to the quantization of vector fields in three-dimensional 
cases. 

For three-dimensional cases, the four components of the 
electromagnetic vector potential are not all independent. So 
we have to fix gauge and to remove those extra degrees of free- 
dom before carrying out the direct quantization of the mass- 
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less vector field. Generally, the quantization of the electro- 
magnetic field in three-dimensional and time-dependent me- 
dia is performed in the Coulomb gauge [il8t 11911 .In these the- 
ories, two unphysical degrees of freedom of the electromag- 
netic field are eliminated from the start and only two trans- 
verse polarizations are quantized. Recently the path integral 
quantization of electromagnetic field for dielectrics in covari- 
ant gauge has been proposed, where all polarizations of the 
photon are quantized [20]. Remarkably, it is found that the 
contributions of the unphysical photons will not be cancelled 
by those of the ghost fields, which implies that there mig ht be 
some physical effects of the unphysical photons ll2ll|22 ^. 

In this work, we try to put forward a Gupta-Bleuler quanti- 
zation procedure to the time-dependent medium. We find that 
in this case, the contributions of the unphysical photons can- 
not be removed completely. Particularly, when the initial state 
is the vacuum state, contributions from unphysical photons 
and transverse photons to the average number are equal. This 
will lead to a detectable surface charge density in the medium 
surface that was not obtained from theories using Coulomb 
gauge ifisl and may be tested experimentally by detect- 
ing the existence of this surface charge density. Moreover, if 
the Gupta-Bleuler theory describes facts more reasonably, it 
will also provide another possible way to detect the dynami- 
cal Casimir effect. 

The article is organized as follows. In Sec. II, we firstly de- 
rive a general Lorentz gauge from Maxwell equations in the 
time-dependent medium. In Sec. Ill, we start by construct- 
ing a new Lagrangian based on this gauge and then present 
a canonical quantization scheme of the electromagnetic field. 
The Gupta-Bleuler condition is introduced in Sec. IV. The di- 
electric permittivity with sinusoidal time dependence is con- 
sidered in Sec. V, and we show that by the rotating-wave ap- 
proximation, it could lead to the effect of unphysical photons, 
namely the surface charge density. Possible implications and 
experimental tests are discussed in Sec. VI. We then conclude 
our work in the last section. 

Throughout this paper, a signature of (+,-,-,-) for the met- 
ric will be adopted. Greek indices are summed over (0,1,2,3) 
while Latin indices denote the spatial components (1,2,3). 
However, letter A is singled out as the polarization state label. 
And natural units c = fi = SQ=fiQ-l will be used. 
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II. GENERAL GAUGE CONDITION 

In this section we will develop a general Lorentz gauge 
adaptive to time-varying media. It is convenient to start with 
the consideration that the pure radiation field and the macro- 
scopic source free Maxwell equations in media take the form 



V ■ D = 0, 



V X E + 5 = 0, 
at 



V ■ B = 0, 



dt 



(1) 



(2) 



(3) 



(4) 



In the following, we shall restrict our discussions to the 
case of linear, nondispersive, time-dependent dielectric me- 
dia, and assume that the constitutive relations between (D, B) 
and (E, H) can be simply expressed by 



gauge condition f^Ao -H V ■ A = o). According to this con- 
straint, Eqs. (|9]l and ( fTOl i decouple into inhomogeneous 
damped wave equations, with the time and space components 
written separately as 

V^Ao - e(f)Ao - 2e(f)^Ao - s(t)—Ao = 0, (12) 

of OP 



V^A- 



di 



52 



0, 



(13) 



D = s{t)E, 



(5) 



where the dot stands for a time derivate. Generally speaking, 
we can see the apparent distinction between these two equa- 
tions, which will lead to different patterns of evolution in time 
and space. In addition, we can attain negative or positive signs 
of the term e(f), which correspond respectively to wave atten- 
uation or amplification. If one makes a special choice e(f) - 1, 
both Eq. ( fT2b and Eq. ( fT3l l will return back to the d'Alembert 
equation. 



III. CANONICAL QUANTIZATION IN 
TIME-DEPENDENT MEDIA 



B = H, 



(6) 



where E(t) is a time-dependent dielectric permittivity. Strictly 
speaking, the dielectric permittivity is complex but in this pa- 
per we assume that it is real ifToll . 

Generally the field strength, E and B, in quantum mechan- 
ics, are regarded as derived quantities and the fundamental 
4- vector A^ = can be introduced according to the re- 

lations 



B = V X A, 



E 



V7. 

-VAo--. 



(7) 



(8) 



With this potential, the two homogeneous Maxwell equations 
(|2]l and ^ are automatically satisfied, and the other two in- 
homogeneous equations ([T]l and ^ can be written in terms of 
Eqs. ^ and ^ as 



(e(f)Ao) ■ 



■V^Ao 



dt\dt 



(e(f)A„) + V ■ A 



d(dA\ 



(9) 



-V -(e(f)Ao) + V-A|. (10) 



From these expressions, we can choose a gauge condition by 
imposing the requirement 



- (e(f)Ao) + V ■ A = 0. 
of 



(11) 



This gauge condition ( fTTT i is a generalization, appropriate to 
the presence of time-dependent dielectrics, of the Lorentz 



In this section we try to canonically quantize the electro- 
magnetic field A'' (x) in time-dependent media. First of aU we 
observe that Eqs. ST% and ( fTlT l can be shown as the equa- 
tions of motion derived from the variational principle with the 
following Lagrangian 

r = -jF.^H''' - ^ {do {s it) A") + diA-f , (14) 



where 



e{t)F''^ F'j 



(15) 



Here the field strength tensor is defined as = d/^Ay - 5vA^ 
and ^ is a parameter that can be chosen freely. To simply 
the calculations, we generally set ^ = 1 which is called the 
Feynman gauge. The extra term on the right hand side of Eq. 
(fl4] i is known as gauge-fixing term constructed by the general 
Lorentz gauge ( fTTT i. We find that this modified Lagrangian 
satisfies canonical quantization prescription. For simplicity 
we may takeA^ - (e(f)Ao,A,). can be shown as the equa- 
tions of motion derived from the variational principle with the 
following Lagrangian (fT4l i as the follows 



d{doA^) 



{-do{e(t)A°),-s(t)doA'). (16) 



For now the II" no longer vanishes and then we can continue 
to carry out the canonical quantization programme. 

As a routine, the quantization of electrodynamics in di- 
electric media is performed by imposing canonical equal-time 
commutation relations (ETCR) 

f A'' (x, f ) , n" (x' , f )] = ig^'d^ (x - x' ) , (17) 
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U'' (x, f) , A" ( x', f)l = [n'' (x, t) , (x', t)] = 0, (18) and 



where g^^ = = diagi l, -1,-1,-1) is the flat Minkowski 
metric. Note that all four components of A'^ and 11'^ have been 
brought into consideration and they obey these relations. 

In order to characterize quantum properties of the field, the 
operator A^ (x) in the Heisenberg picture is expanded in the 
polarization vectors e^"^^ as 



^(^)t 



(19) 

where a^'*' (f) and a*'^" (f) are the annihilation and creation op- 
erators. For each k, the four polarization 4-vectors ef\ e'P, 

efj^ and efj^ are linearly independent, and are usually chosen 
to be real. Here (ph (x) and ujtt are respectively instantaneous 
mode functions and eigenfrequencies, which obey the equa- 
tion 



(v2+£(f)w2,)0„(x) = 0. 



(20) 



The index f is regarded as a parameter indicating that these 
mode functions and eigenfrequencies depend on the value of 
s (f) at each time f . 

In what follow we will confine our consideration to the elec- 
tromagnetic field inside a three-dimensional cavity filled with 
time-dependent media. By choosing an orthonormal set of 
mode functions {^^, (x)}, they are required to satisfy the or- 
thonormal conditions lfl2ll 



1 £xB (f) 4>kt (X) 4>l,i (X) = 6kk' , 

Jv 



(21) 



where the integration is performed over volume V. In obtain- 
ing Eq. (|5]l we have assumed that the medium is isotropic so 
that the mode function may be described by plane wave 



eigenfrequencies 



can be expressed as wj;, = 



At present, to exhibit the obvious physical meanings of the 
expansion (fT9] l, the polarization vectors el^' can be chosen in 



a fixed reference frame as |2 



( 1 ^ 




( ^ 




i ° 1 







ki 


1 (2) 


k\k-i 







-ki 


— - 

h •'' 




loj 




{ J 




\ "ii ) 



1 

kk\\ 



k,fi 




ki 
ki 
h ) 
(26) 



1 



with k\\ = yjk^ + k\ and k - ^k^^ + ky These vectors form a 

set of basis vectors that satisfy the following orthonormal and 
complete relations 



AA' 



(-1) W) _ 



(27) 



where g'^'^ - diag (1,-1,-1,-1), like g^^^ . It is seen from Eqs. 
and (l26b that the polarization vectors with A - 1 , 2 are 

while those with /I = 0, 3 



(0) 



-e-e 



(3) 



transversal and satisfy { ■ e'^ 
are timelike and longitudinal constrained by ^ ■ e 

Substituting Eqs. ( |26] | and ( |22] | into the expansion ( fT9] l we 
obtain the standard form 



1 . . ^-ik.x^ 



+eflaf^ (t) e-i (f) + (f) g-J (f) 

+eflafUt)E-Ht)} + c.c., (28) 

with j - 1,2,3. Here the polarizations A = 1,2 correspond 
to the transversal parts of the potential (l2ll which can be also 
given in Coulomb gauge. We then insert Eq.(l28ll into the wave 
equations (fT2] l and (fT3T i. and obtain the evolution equations of 

a[^'> (f) and af'^' (f) 



(pkt (x) = L 



-3/2 



1 



ikx 



(22) 



with the length of one side of the cavity L. Furthermore, the 
periodic boundary conditions now are read as 

27m_, 2nny 27m, 

k\ — , kj — -, k-i — , n,., «,„ n, = U, ±1, +2, ■ ■ ■. 

L L L " ■" - 

(23) 

Note that in the formal limit L — > oo, these wave vectors 
become continuous rather than discrete, which describes the 
electromagnetic field in free space filled with media. Then we 
can define a 4-vector 



(Ve (0^0, ^1,^2, ^3), 



(24) 



with kd - iOkf From dSOll, (l22]i and (|24li it follows that the 
energy-momentum relation of photons can be simply written 
as 



(25) 



daf (f) 6-""^'' 
dt 

da[^^Ht)e"'J'' 
dt 



-itOkia^ {t)e ' +- a, (t)e ' 

Z OJki 



(29a) 



2ojk, 



(29b) 



where the upper sign (-n) in the right hand side of the equations 
denotes timelike photon A = 0, and the lower sign (-) denotes 
transverse and longitudinal photons A = 1,2,3. Now, at the 
initial time ? = 0, we assume a^^ (0) = a^'*^ to be usual annihi- 
lation operators. If we set the parameter e (f) - 1, then all co- 
efficients dtitt are zero. In this special case, Eq. (I29al i reduces 

to —^7^ - -icJktd^k ^ ({)■ Hence the time evolution of the an- 
nihilation operators is simply expressed by a^'*' (t) = a^f'e^"^''' 
with a time-independent frequency a>k- Similarly, we can get 
the time evolution of the creation operators in the same way. 
That is to say, they have been reduced to the expression in 
vacuum |23| . 
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Now we put Eq. ( l28T l into the commutation rules (fTTI i and 
( fTST l. and then derive the commutation relations for the com- 

' "(f) as 



ponents a^"^^ (f) and a^''^ 



I' 



,(-1) 



(f),«r'(o] 



J-''*' 5^ (k ■ 



(30) 



[fl^-*' (f) , fl^f > (f)] = [fll'^"'" (f) , fllf ' (f)] =0. (31) 

It can be seen that for longitudinal and transverse photons 
A = 1,2,3 these relations are subject to the ordinary bosonic 
quantization scheme. Nevertheless, for timelike photons A - 
there is an extra minus sign on the right-hand side, which 
would bring about a remarkable consequence that the Hilbert 
space of the photon field carries an indefinite metric. 

Then, for the operators aj,'' (0 and a^'*^' (f), Heisenberg 
equations of motion are given by 

iafHt) = [afUt),H,ff]. (32) 

Comparing the above equations ( |32] | with the evolution equa- 
tions ( |29] |. the effective Hamiltonian operator in the Heisen- 
berg picture is found to be 



IV. THE GUPTA-BLEULER METHOD IN 
TIME-DEPENDENT MEDIA 

So far we have not been really dealing with Maxwell's the- 
ory in media, since we have adjusted the Lagrangian (fT4l i. To 
recover it, let us now turn our attention to the general Lorentz 
gauge ( fTTI ). However, after the quantization we cannot simply 
take this gauge condition as an operator identity. Comparing 
with the classical Gupta-Bleuler method ll27[|28ll . we impose 
the general gauge condition on physical states \(p) as the ex- 
pectation value 



[do{s(t)A^) + djAj]\<p}^0. 



(37) 



For simplicity, the above formulas can be equivalently ex- 
pressed as 



{do{s{t)A'') + djAjf^ 10) =0, 



(38) 



where only the destruction part of this gauge condition is al- 
lowed to act on physical states. Now, we insert ( |28] l into ( |38] l 



and take into account of { ■ e. 



(3) 



4;'=''2^ = 0andref = 
Then, the constraint at the initial time can be rewritten as 



k ,1=0 



i k (t) 



(33) 



with K (t) - Ve(f). This Hamiltonian in the Schrodinger pic- 
ture is given simply as ll24ll 



k "k 



i k it) / 

'4K(t) 



— a 



f) 



(34) 



We note that in the formal Umit e(f) — > 1 the Hamiltonian 
( l34l l turns into the Hamiltonian for the electromagnetic field 
in vacuum |23]. For performing the perturbation theory we 
would like to further investigate the system's evolution in 
the interaction picture. Through the unitary transformation 
t/o (0 = exp (-/ ^' c/t/Zq (t) , we get the system Hamiltonian 
in the interaction picture as j 25l 

(35) 
and 



where the free part is Hq 
interaction Hamiltonian is 



Ijk 2ji=o (-^,1,1) Wfefl^ fl^ 



Hi 



ZjZj4^(f) V 



),(36) 



k i=0 



with Qi (f) - cokrclT. It is seen that the first term in Eq 

is two-photon process characterized by the term a^'*'^ fl[^' ' , so 
that photon pairs can be created from the vacuum state. Sec- 
ondly, in contrast to f25','2^, we are interested here in whether 
it is likely to create longitudinal and timelike photons, in ad- 
dition to transverse photons. This will be discussed in detail 
in Sec.V. 



0, 



(39) 



where |0)o is the initial quantum state (at t — 0). This ex- 
pression provides a link between the timelike and longitudinal 
photons. And from Eq. ( |39] l and its adjoint we have 



0<<^|flf ^ flfV)o - 0< 



= 0. 



(40) 



By this means we will remove the contributions of the unphys- 
ical photons to the expectation value of the number operator 
and energy in the initial time. However, with the time evolu- 
tion, we will find that in the medium the contributions of the 
longitudinal and timelike photons cannot be cancelled out by 
each other In this case, it is observed that there is a surface 
charge density in the medium surface. We will devote a full 
discussion to this subject in the next section. 



V. SURFACE CHARGE DENSITY 

In the following, by means of dielectric permittivity chang- 
ing in periodic time dependence, we will study the nature of 
electromagnetic fields in the medium and investigate the pos- 
sible impact on the detection of the dynamical Casimir effect 
. Now we assume sinusoidal time-varying permittivity 



e{f)^e + 26 sin {2Q.t), 



(41) 



where e is a constant, the amplitude 6/e <si 1, and Q. - a>i 
is the static value of the fundamental mode for e (f ) = £ fl^ . 
Then by taking dielectric permittivity (l4]linto the Hamilto- 
nian operator i/J ( [36] l and using the rotating-wave approxima- 
tion [25. .261 . we can obtain 



H 



4e 



,1=0 



)- 



(42) 
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Generally assuming that the electromagnetic field is initially 
in thermal equilibrium states at finite temperature T ll26ll . the 
expectation value of the number operator in the A;th mode can 
be written in the form 



{Nk (0> * i-gAA) { exp [iH[t)af'af exp (-///{f))o 
+5Hsinh2(^Qf)<flf'-4^'+4«flf')o) 



,)=0 



^=0 



(43) 



Because of the formula gOli, <A^f *)o + (A^f >o = 0, we get 



2 2 „ 



2sinh^(|or) 



^=0,3 



(44) 



The first term Ya=\{N^]^^)Q denotes the initial photons contain 
only the transverse photons, while the second term represents 
the average number of photons created grows exponentially 
with time. It is noticeable that here the contributions of the 
longitudinal and timelike photons are nonzero. Nevertheless, 
only the transverse part of photon generation A - 1,2 has 
been enhanced by the initial thermal photons. This means that 
the number of created unphysical photons /I = 0, 3 is much 
smaller than that of the transverse photons at nonzero temper- 
ature. 

As a next step, we let the vacuum state |0) be the initial 
quantum state (at T - 0), then the total number of photons 
created from vacuum in the time f 



A—0 



(45) 



According to this result, we find that all polarization compo- 
nents have the same contribution to the number of photons 
generated from vacuum. This is to say, the number of un- 
physical photons is equal to those of transverse photons in 
this case. 

Further, we would exhibit the content of the electric and 



magnetic fields. Through £, (x) - 



jA; - we can obtain 



from Eq. 
Ei(x)-. 




y/2k 



/(-«f(0 + af (0) 



(1) / X 

-aY> (f) + 



( hh ^ 
-kl 



(46) 



and from B, (x) = eij^djAk for the magnetic field 

ik 



Bi(x) 



( kiki ^ 
kiki 
-kl 



-kl 




—a?' (t)} + c.c. 



(47) 



It is seen that longitudinal and timelike photons exist only 
in the longitudinal component of the electric field, not in the 
transverse electric and the magnetic fields. 

Therefore we only need to investigate the effects of the lon- 
gitudinal electric field. In contrast to ll29ll . the energy of the 
unphysical photons can be represented from (l4&t as 

i r d'xe (f) Ej, = Yj ^-tr [af (0 af (f) - flf (f) (t)] , 
k 

(48) 

where £// is the longitudinal part of the electric field operator 
E; (x). At the same time, from [30] this longitudinal electric 
field on either side of the boundary surface satisfies the bound- 
ary condition 



(«(Oe^;)-e;; 



'acmim 

II 



(49) 



where cr is the surface charge density. In vacuum, the lon- 
gitudinal component of the elective field E™'"'""' is null 
Finally, using Eq. (|45] | and the expectation value of (|48]) . we 
get 



cr =s 2 sinh 



.1/4 



L2 



(50) 



Note that when the dielectric permittivity varies with time, 
the contributions of the longitudinal and timelike photons will 
no longer vanish, which lead to a detectable surface charge 
density. 

This result has not been obtained from the previous quan- 
tum theory ifisi [T9I1 . Thus, in time-dependent media, the 
Gupta-Bleuler quantization scheme and the quantization in 
Coulomb gauge are no equivalent. Theoretically it is impos- 
sible to tell which quantum theory is more reasonable, but ex- 
perimentally this may be tested through the measurement of 
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the surface charge density. If the Gupta-Bleuler theory is rea- 
sonably in this case, we can also propose another experimental 
method for the detection of the dynamical Casimir effect. 

In the experiment fis", 16"], one may accomplish the wall 
motion by manipulating the dielectric permittivity of the 
medium in the cavity. And from above mentioned reviews 
the unphysical photon created from time dependent media 
may lead to a surface charge density. Let us insert some ex- 
plicit number, and estimate this surface charge density. For 
instance, we use the same experimental parameters consid- 
ered in Ref. lfT6ll : 6/e ~ 10"^ and the frequency of the low- 
est mode is wi ~ GHz for Lq ~ O.lm, then approximately 
cr X lO^'^C ■ m^^ will be created in the cavity during 1 s. 

Though detecting this surface charge density, we can verify 
whether the unphysical photons make a nonzero contribution 
to observable quantities in time-dependent media. Namely, 
are the quantization in the Lorentz gauge and that in the 
Coulomb gauge completely equivalent in this case? 

At the same time, from the DCE 1 16] detection experiment, 
it is difficult to distinguish between the background of thermal 
photons and the creation of transverse photons. Therefore, we 
suggest to detect the surface charge density, hopefully it will 
provide another probability to observe the dynamical Casimir 
effect in the near future. 



VI. CONCLUSIONS 

Starting from the source free Maxwell equations in media, 
we presented a gauge condition, which is a generalization of 
the Lorentz gauge condition appropriate to time varying me- 
dia. Two different wave equations were obtained in the direc- 
tions of time and space. It is important to notice that e(/) is in 



general nonzero, which may lead to photon generation even 
from vacuum. 

Through the general gauge (fTTT i a modified Lagrangian was 
given in Sec. Ill, which makes II" no longer vanish, and car- 
ries out the canonical quantization procedure. Assuming that 
the mode functions can be described by plane wave, and using 
the polarization vectors e^"^^, the standard expansion of the op- 
erator (x) and the evolution equations of a^'^' (f) and a^'*^' (f) 
were deduced. Based on the Heisenberg equations of motion, 
the effective Hamiltonian was immediately derived. But this 
Hamiltonian formula ( |33T l suggests that all polarization com- 
ponents contribute to the photon creation. In accordance with 
the Gupta-Bleuler method in time-dependent media, imposing 
the constraint condition ( fTTT i on the physical state, we were 
able to get rid of contributions from the longitudinal and the 
timelike photons at the initial time. 

Finally, by assuming sinusoidal time varying permittivity, 
we were able to simulate this effect. It is shown that not only 
the transverse photons but also the longitudinal and timelike 
photons were obtained in the experiments of DCE. Particu- 
larly, for an initial vacuum state, all polarization components 
have the same contribution to the number of photon genera- 
tion. Therefore, we can get a detectable surface charge den- 
sity in the medium surface, and provide another experimental 
trick for confirming one of the fundamental predictions of the 
DCE. 
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